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In order to overcome the limitations of small system sizes in spin-glass simulations, we 
investigate the one-dimensional Ising spin chain with power-law interactions. The model 
has the advantage over traditional higher-dimensional Hamiltonians in that a large range of 
system sizes can be studied. In addition, the universality class of the model can be changed 
by tuning the power law exponent, thus allowing us to scan from the mean-field to long-range 
and short-range universality classes. We illustrate the advantages of this model by studying 
the nature of the spin glass state where our results hint towards a replica symmetry breaking 
scenario. We also compute ground-state energy distributions and show that mean-field and 
non- mean-field models are intrinsically different. 

§1. Introduction 

Because spin glasses generally equilibrate very slowly due to their glassy behav- 
ior, simulations can usually only be performed on small system sizes. The intro- 
duction of novel algorithms, such as exchange Monte Carlo (parallel tempering) 
have helped considerably to overcome these system-size limitations. Nevertheless, 
in general, only modest system sizes can be simulated at temperatures well below 
the critical temperature. In what follows we introduce the one-dimensional Ising 
chain with random power-law interactions in order to simulate large length scales. 
We illustrate the advantages of this model with two applications: the nature of the 
spin-glass state and ground-state energy distributions in spin glasses. 

§2. Model & Numerical Method 

The Hamiltonian for the one-dimensional (ID) long-range Ising spin glass with 
power-law interactions is given by*"™*" 



H = -J2 J ijS i S j , Jij-^r, 7>(e ij ) = -±=exp(-<? j /2), (2-1) 

i,3 13 



where the sites % lie on a ring of length L to ensure periodic boundary conditions, 
and Si = ±1 represent Ising spins. The sum is over all spins on the chain and 
the couplings Jy are Gaussian distributed, and divided by the geometric distance 
between the spins to a power a. The model has a rich phase diagram in the d—a 
planed 1 For a < 0.5 the model is in the mean-field phase with an SK universality 
class, for 0.5 < a < 1.0 the model has long-range critical exponents with a finite 
transition temperature T c , whereas T c = for 1.0 < a < 1.75. For a > 1.75, again 
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T c = but in a short-range universality class. In addition, there is a prediction 
from droplet argument JSJ'EJ that 9 = d — a, where 9 is the stiffness exponent for 
domain- wall excitations. 

For all of our simulations, we use the parallel tempering Monte Carlo methocP'^ 
as it allows us to study larger systems (up to L = 512) at very low temperatures. 
Details about the simulations as well as equilibration tests can be found in Ref . EJ) JHI . 

§3. Applications 

3.1. Nature of the spin-glass state 

Two main theories attempt to describe the nature of the spin-glass state: The 
replica symmetry-breaking (RSB) picture and the "droplet picture" (DP). RSB pre- 
dicts that droplet excitations involving a finite fraction of the spins cost only a finite 
energy in the thermodynamic limit, and that the surface of these excitations has a 
fractal dimension d s equal to the space dimension d. In the DP, excitations have an 
energy proportional to £ s , where i is the characteristic length scale of the droplet 
and 9 is a positive stiffness exponent. In addition, the excitation surfaces are frac- 
tal with d s < d.*' . Differences between the various pictures can be quantified by 
studying^'^ P(q), the distribution of the spin-glass overlap q = L~^Y^=i^f^, 
where "a" and refer to two replicas of the system with the same disorder. The 
RSB picture predicts a nontrivial distribution with a finite weight in the tail around 
q = 0, independent of system size. In contrast, the droplet picture predicts that 
P(q) is trivial in the thermodynamic limit with only two peaks at ±(?EA) where q^A 
is the Edwards-Anderson order parameter. For finite systems there is also a tail 
down to q = 0, which vanishes in the thermodynamic limit like*'»^' P(0) ~ L~ e 
with & = 9**1 

Figure n left panel, shows data for P(q) at a = 0.75 and T = 0.10, well be- 
low T c PS 0.62.ENn) There is a peak for large q and a tail down to q = that is 
independent of system size. A more precise determination of the size dependence 
of P(0) is shown in the inset of Fig. ^ left panel, where, to improve statistics, we 
average over q values with \q\ < 0.50. The expected behavior in the droplet model is 
P(0) ~ L~ 6 ', with 9' = 9 wherPS 9 = d-a. The dashed line has slope -0.25, the 
expected value for a = 0.75 according to the droplet model. The size dependence 
is consistent with a constant P(0), which implies that the energy to create a large 
excitation does not increase with size, at least for the range of sizes studied here. In 
Ref. EI) we estimate the fractal dimension d s of the system-size excitation surfaces 
and we find that d s > 0.95 ± 0.05, i.e., d s « d, in agreement with RSB. 

Since there are analytic predictions for the stiffness exponent of the ID chain 
(9 = d — a), and the droplet picture predicts 9 = 9', we compute 9 directly from 



*' Recently Krzakala and Martin, as well as Palassini and Young, suggest an intermediate picture 
in which droplets have a fractal surface, and their energy is finite in the thermodynamic limit. See 
Ref. [SJ and references therein for details. 

**' Note that we explicitly distinguish between 9' , the stiffness exponent for droplet excitations, 
and 8, the exponent for domain walls, because DP predicts 9 = 9', which we want to test here. 
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Fig. 1. Left panel: Distribution of the spin overlap P(q) at T = 0.10 and a — 0.75 (T c finite, 
long-range universality class) for several system sizes L. The data are independent of L at 
g ~ 0, indicating that 6' = 0. This is shown in detail in the inset where P(0) ~ L~ e is plotted 
as a function of L together with the prediction from the droplet model (dashed line with slope 
—0 = —(1 — a) = —0.25). Right panel: stiffness exponent 6 as a function of a at T — from 
domain-wall calculations. The data follow well the prediction of the droplet picture 6 — d — a 
(solid line). The star ("*") marks a — 0.75, where 6 = 0.173 ± 0.005, which disagrees with 
0' — from finite-temperature calculations (left panel). The dashed lines are guides to the eye. 
Note also that 6 — > 0.3 for a — > 0, the Sherrington-Kirkpatrick (SK) model limit. 



zero-temperature domain-wall calculations for different values of a using parallel 
tempering as an optimization algorithm.^ EJCU The change in energy AE induced 
by a change in boundary conditions from periodic (P) to antiperiodic (AP) scales as 
AE = [\E AP - E P \ } av ~ LP , where [• • -] av represents a disorder average. Data for 
0{o~) are shown in Fig. ^ right panel. The solid line represents the prediction from 
the droplet model, well followed by the data. For a = 0.75 ("*" in Fig. ^) there is a 
clear positive stiffness exponent: 9{a = 0.75) = 0.173 ± 0.005. This is in contrast to 
9' = from finite-temperature simulations (Fig. left panel) . Therefore we find a 
disagreement with the droplet model in that 9' ^ 9 for o large range of sizes. 

3.2. Ground-state energy distributions 

There has been recent interest in how ground-state energy distributions scale in 
spin glasses. Considerable work has been done for the Sherrington-Kirkpatrick (SK) 
modeP^ 1 where evidence for skewed energy distributions, well fitted by a modified 
Gumbel distribution^ with m = 6, is found. In order to test if these results are 
intrinsic to the mean-field SK model we compute ground-state energy distributions 
for the ID chain for different system sizes and different values of a. The advantage 
of the ID chain is that the crossover from mean-field to short-range behavior can be 
probed for a large range of sizes. 

In Fig. left panel, we show rescaled ground-state energy distributions P{E) 
with mean (E) and standard deviation oe- The data show a clear asymmetry. In 
order to quantify these effects we show data for the skewness C,e of the distributions 
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Fig. 2. Left panel: Rescaled ground-state energy distributions for a — 0.00 (SK limit) as a function 
of system size L (10 samples). The data are clearly asymmetric; the dashed line is a guide to 
the eye. Right panel: Skewness £g as a function of L for different values of a. In the mean-field 
region (a < 0.5) the C,e — > const, for L —> oo. For a > 0.5 the distributions become symmetric 
in the thermodynamic limit (£e ~ L 1 , 7 < 0). 



as a function of system size L for different values of the power-law exponent a in 
Fig. |2 right panel. One can see that for a < 0.5, where the model exhibits mean-field 
behavior, the skewness tends to a constant in the thermodynamic limit, indicating 
that the skewness of the distribution persists even for infinitely large system sizes. 
For a > 0.5, the skewness decays with a power law of the system size, indicating that 
in the non-mean field region the ground-state energy distributions become symmetric 
in the thermodynamic limit. This shows that mean-field and non-mean field models 
are expected to behave differently, and that, in particular, intrinsic length scales in 
the mean-field model seem to scale with system size. 



§4. Conclusions 



By using a one-dimensional Ising spin glass with power law interactions, we 
have been able to overcome the usual limitation of small system sizes in spin-glass 
simulations. We illustrate the advantages of the model on the nature of the spin- 
glass state where we show that droplet and domain-wall excitations scale differently 
for large system sizes, indicating that the droplet model cannot be correct, at least 
for the ID chain. In addition, we compute the ground-state energy distributions 
for the one-dimensional Ising chain and show that, in the mean-field regime, energy 
distributions remain skewed in the thermodynamic limit indicating that mean-field 
and non-mean-field models behave differently. Our results are in agreement with 
previous results for the SK model (cr = 0). In the future, we intend to study energy 
distributions at finite temperature, as well as the effects of degenerate ground states. 
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